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Abstract 

Alfvenic waves are thought to play an important role 
in coronal heating and solar wind acceleration. Here we 
investigate the dissipation of such waves due to phase 
mixing at the presence of shear flow and field in the 
stratified atmosphere of solar spicules. The initial flow 
is assumed to be directed along spicule axis and to vary 
linearly in the x direction and the equilibrium mag- 
netic field is taken 2-dimensional and divergence-free. 
It is determined that the shear flow and field can fas- 
ten the damping of standing Alfven waves. In spite 
of propagating Alfven waves, standing Alfven waves in 
Solar spicules dissipate in a few periods. As height in- 
creases, the perturbed velocity amplitude does increase 
in contrast to the behavior of perturbed magnetic field. 
Moreover, it should be emphasized that the stratifica- 
tion due to gravity, shear flow and field are the facts 
that should be considered in MHD models in spicules. 

Keywords Sun: spicules • Alfven waves: phase mixing 
• shear flow • shear field 



1 Introduction 

Phase mixing has been proposed as a mechanism of ef- 
ficiently dissipating Alfven waves in the solar corona by 
Hevvaerts fe Priest) jl983l) . iKarami fe Ebrahimil |2009l) 
calculated numerically the damping times of standing 
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Alfven waves in the presence of viscosity and resis- 
tivity in coronal loops. They con cluded that the ex- 
ponen tial damping law obtained by iHevvaerts fc Priest 
( 19831) in ti me is valid for the L undq u ist nu mbers higher 
than 10 7 . be Moortel et al.l |l999l l2000l ) studied the 
effect of stratification and diverging background mag- 
netic field on phase mixing, and found that the wave- 
length of an Alfven wave is shortened as it propagates 
outwards which enhances the generation of gradients. 
They concluded that the convection of wave energy 
into heating the plasma occurs at lower heights than 
in a uniform model. Moreover, the combined effect of 
stratification and diverging background m agnetic field 



depen ds on the geometry of configuration. ISmith et al 



(|200?T) showed that the enhanced phase mixing mech- 
anism can dissipate Alfven waves at heights less than 
half. Moreover, it can occur in divergent and stratified 
coronal structures only when the ratio of magnetic and 
density scale heigh ts is lower than 2. 
Kaghashvili ( 19991 ) studied the effect of inhomogeneous 



flow on converting the Alfven waves into other types of 
MHD waves that can dissipated efficiently. It is found 
that in the divergent geometry of magnetic fields, this 
mechanism becomes essential in the rapidly expand- 
ing re gions where t h e velo city shear flow is relatively 



large. ISaleem et al.l ([20071 ) studied drift modes driven 



by shear plasma flow in spicules and concluded that 
the density inhomogeneity and the shear in the plasma 
flow should be expected as facts in their atmosphere. 
These waves propagating on spicules can be efficiently 
dissipated in the regions with sheared magnetic field in 
such a way that can heat the solar corona. 
Spicules appear as grass-like, thin and elongated 

structures in images of the solar lower a tmosphere 

Zaaarashvili fc ErdelvJ2009[ISterlingll200(i . iKukhianidze et al.| 



20061) : IZaaarashvili et all ( 2007 ) by analyzing the 



height series of Ha spectra concluded that transverse 
oscillations can be caused by propagating kink waves in 



spicul es. De Pontieu et al. ( 2007 ); Okamoto k, De Pontieufcfily ionized H plasma is 2.2 x 10 17 T 5 / 2 kg m 1 s 1 



( 20111 ) by making use of Hinode observations inter- 
preted that transverse oscillations of spicule axis are sig- 
natures of Alfven waves, and concluded that standing 
Alfven waves occurred at the rate of about 20%. More 



recently lEbadi et al.l (|2012al ) based on Hinode /SOT 
observations concluded that spicule axis oscillations 
are closed to the standing pattern as the y do not see 
any up ward or downward propagation. lEbadi et al 



d2012bl) Dcrformed the phase mixing model of propa- 
gating Alfven waves in spicule conditions and concluded 
that the damping times are much longer than spicule 
lifetimes. These results motivated us to study the phase 
mixing of standing Alfven waves in a stratified atmo- 
sphere in the presence of shear flow and field. To do 
so, the section 2 gives the basic equations and the the- 
oretical model. In section 3 the numerical results are 
presented and discussed, and a brief summary is fol- 
lowed in section 4. 



2 Theoretical modeling 

In the present work we keep the effects of stratification 
due to gravity in 2D x-z plane in the presence of shear 
flow and shear field. The phase mixing and dissipation 
of standing Alfven waves in a region with non-uniform 
Alfven velocity are studied. MHD equations, governing 
the plasma dynamics are as follows: 



dp 
dt 



V • (pv) = 0, 



(1) 



/°^-+p( v 'V)v = — Vp + pg+ — (V x B) x B + pi/V 2 v, 
at (x 

(2) 



— - = V x (v x B) + ??V 2 B, 
at 



pRT 



(3) 



(4) 



( Priestl 119821 ). We assume that the spicules are highly 
dynamic with speeds that are significant fractions of 
the Alfven speed. The perturbations are assumed in- 
dependent of y, with a polarization in y direction, i.e.: 



B 



vo{x)k + v y (x, z, t)j 
B 0x (x, z)i + B 0z (x, z)k 



■b v (x,z,t)j. 



(6) 



The backgroun d flow is assumed t o vary linearly in 
the x-direction as (Poedts et al. 19981): 



v = v (x) k = (x- l)v k, 



(7) 



and the equilibrium sheared mag netic filed is two- 
dimensional and divergence-fre e as (jSmith et al.l 12007 
Murawski fc Zaaarashvili 2010l) : 



B 0x (x,z) 
B 0z (x,z) 



= B e- kBZ sm[k B {x - 
= B e~ kBZ cos[k B (x 



1)] 
■1)], 



(8) 



Therefore, the pressure gradient is balanced by the 
gravity force, which is assumed to be g=-g k via this 
equation: 



-Vp + Pog = 0, 
and the pressure in an equilibrium state is: 



Pa = Po(x)e 



-z/H 



(9) 



(10) 



Moreover, the equilibrium density profile is given the 
form: 



po = p (x)e z/H , 
with 



H = 



RT 

pg 1 



(ii) 



(12) 



where H is the pressure scale height. Taking into ac- 
count these assumptions, the linearized dimcnsionless 
form of Eqs. [U gj and H yield: 



V • B = 0, 



(5) 



where v and r\ arc constant viscosity and resistivity co- 
efficients, and other quantities have their usual mean- 
ings. In particular, typical values for r\ in the solar 
chromosphere and corona are 8 x 10 8 T -3 / 2 m 2 s -1 and 
10 9 T -3 / 2 m 2 s _1 , respectively. The value of pv for a 



and 



i x db v „ . . db v 
B x(x, 2)—^ + B 0z (x, z)— It- 
ox oz 



( \ dv v 
oz 



(13) 
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Ok - U ( x zA+Bn (x z)*± 
-vo(x)^- + i]V 2 b y , 



(14) 



where velocities, magnetic field, time and space co- 
ordinates are normalized to Vao = Bo / \Jppo (with 
po as the plasma density at z = 0), Bo, r (Alfvcn 
time), a (spicule radius), respectively. Also the re- 
sistivity and viscosity coefficients arc normalized to 
a 2 jr. The first and second terms in the right hand side 
of Eqs. [T3] and [14] representing effects of shear fields 
and flows, respectively. Va(x,z) is the Alfven veloc- 
ity, which for a phase mixed a nd stratified atmosphere 
due to gravity is assumed t o be (|De Moortel et al ll999t 
Karami fc Ebrahimil [20091) : 



V A (x,z) = V A0 e z/2H [2 + t'Anh[a(x - 1)]]. 



(15) 



Here, parameter a controls the size of inhomogeneity 
across the magnetic field. 

The set of Eqs.[l3l and [T4l should be solved under these 
initial conditions: 



v y (x, z, t = 0) 



Vao exp 



1 ,X — l.n 
(" 

2 



b y (x,z,t = 0) = Asm.(nx) sin(7rz), 



sm(kz)e z/4H 
(16) 



where d is the width of the initial packet and A = 10 -7 . 
Figure [T] is the plot of initial wave packet and magnetic 
field given by Eq. [16] for d = 0.3a (a is the spicule 
radius). Also, the boundary conditions are taken to be: 



v y (x = 0, z, t) = v y (x = 2, z, t) = 
b y (x = 0, z, i) = b y (x = 2, z, t) = 0, 

and 



v y (x, z = 0, t) = v y (x, z = 8,t) 
b y (x, z = 0, t) = b y (x, z = 8, t) - 



0. 



(17) 



(18) 



Boundary conditions in Eq. [TS] imply that we are 
dealing with standing waves. 



v (x,z,t=0) 




b (x,z,t=0) 




Fig. 1 (Color online) The initial wave packet and mag- 
netic field in x — z space are presented from top to bottom, 
respectively. 



3 Numerical results and discussion 

To solve the coupled Eqs. [T31 and [TU numerically, the 
finite difference and the Fourth- Order Runge-Kutta 
methods are used to take the space and time deriva- 
tives, respectively. We set the number of mesh grid 
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10 11 CTO 3 , 

^/c" 1 (uni- 
= pi/ 8 (di- 



points as 256 x 256. In addition, the time step is cho- 
sen as 0.001, and the system length in the x and z 
dimensions (simulation box sizes) are set to be (0,2) 
and (0,8). The parameters in spicule env i ronment are 
as foll ows ( Murawski fc Zaaarashvililhoiot Ebadi et"all 
l2012al) : a (spicule radius)=1000 km, d = 0.3a = 300fcm 
(the width of gaussian paket), L=8000 km (Spicule 
length), t>o = 25km/s, Bq = 10G, n e = 
T = 8000.fi:, g = 272ms" 2 , R = 8300to 2 s 
versal gas constant), Vao = 50km/ s, k 
mcnsionlcss wavenumber normalized to a), ks = pi/16, 
v = 10 3 m 2 s" 1 , T] = 10 3 m 2 s- 1 , p, = 0.6, H = 500 km, 
t = 20 s, and a = 2 (jOkamoto fc De Pontieull201ll ). 
Figure [5] shows the perturbed velocity variations with 
respect to time for x = 1000 km, z = 1300 km; x = 1000 
km, z = 4000 km; and x = 1000 km, z = 6700 km re- 
spectively. We presented the perturbed magnetic field 
variations, obtained from our numerical analysis, in 
Figure [3] for x = 1000 km, z = 1300 km; x = 1000 
km, z = 4000 km; and x = 1000 km, z = 6700 km 
respectively. In these plots the perturbed velocity and 
magnetic field arc normalized to Vao and Bq respec- 
tively. In each set of plots it is appeared that both the 
perturbed velocity and magnetic field are damped at 
the middle stage of phase mixing. 

As height increases, the perturbed velocity ampli- 
tude does increase in contrast to the behavior of per- 
turbed magnetic field. This means that with an in- 
crease in height, amplitude of velocity oscillations is 
expanded due to significant decrease in density, which 
acts as inertia against oscillations. Similar results are 
observed by time-distance an alysis of Solar spicule os- 
cillations (jEbadi et al.ll2012al) . It is worth to note that 
the density stratification influence on the magnetic field 
is negligible, which is in agreement with Solar Opti 



cal T elescope observations of Solar spicules (jVerth et al 



2011). 



Figures Hl[5] show the contour plots of the perturbed 
velocity and magnetic field with respect to x, z for 
t = lOr, t = 30t, and t = lOOr. They show that in 
the presence of stratification due to gravity, and shear 
flow and field the damping takes place in time as an im- 
portant point of these graphs. In other words, in spite 
of propagating waves, standing waves arc dissipated af- 
ter a few periods due to phase mixing in the presence 
of shear flow and field. Once again, increment (decre- 
ment) of the perturbed velocity (magnetic) amplitude 
with height at different times are obvious from these 
figures. On top of these, latter graphs show stochas- 
tic evolution pattern regarding the x-direction. In fact, 
taking a package of Alfvenic waves with different speeds 
at different 'x' points results in such non-uniform dis- 
tribution of oscillations. To be more rigorous, consider 
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Fig. 2 The perturbed velocity variations with respect to 
time in x = 1000 km, z = 1300 km; x = 1000 km, z = 4000 
km; and x — 1000 km, z = 6700 km respectively from top 
to bottom are showed. 
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Fig. 3 The same as in Figure [2] but for the perturbed 
magnetic field. 



Eqs.[l3l and[T4]where the z-derivatives are omitted and 
main terms are kept only. For instance, one gets: 



dt 2 



Vl(x,z) 



Box(x, z) 



dB Qx dv y 
dx dx 



B 



2 &\ 

0x dx 2 



(19) 



for the perturbed velocity field. Differential Equa- 
tion \T9\ simply, implies that v y depends to 'x' but un- 
fortunately not in the form of simple function. Addi- 
tionally, having ignored the z-derivatives, still, v y varies 
with V through Va(x, z) and Bq x (x, z) terms. Note 
that, although, one can discuss about the above point 
in detail, but, here, we are more interested in study- 
ing the effects of phase mixing on time evolution and 
propagation of Alfvenic waves in the propagation di- 
rection rather than looking for variations in the specific 
perpendicular 'x' direction. 

To estimate the damping time it is suitable to cal- 
culate the total energy (kinetic energy plus magnetic 
energy) per unit of length in y direction as: 



E t ot{t) 



16tt 
B 2 aL 



dx I dz[p(x, z)v y (x, z) + b y (x, z)]. 

(20) 



In Figure [5] the kinetic energy, magnetic energy, and 
total energy normalized to the initial total energy are 
presented respectively from top to bottom. It is clear 
from energy plots that the standing Alfven waves un- 
der spicule conditions can dissipated after a few pe- 
riods. This is in agreement with the fact that spicules 
have short lifetimes, and are transient phenomena. The 
energy flux, stored in oscillating spicule axis, is of the 
order of coronal energy loss in quiet sun. Therefore, dis- 
sipation of standing Alfven waves can count as a candi- 
date mechanism in transferring and as a result heating 
the solar corona. Moreover, it should be emphasized 
that the stratification due to gravity, shear flow and 
field are the facts that should exist in MHD models in 
spicules. 



4 Conclusion 

In our model, we assume that spicules are small scale 
structures with an initial shear flow and field, and a 
uniform temperature along them. Density variation 
along the spicule axis is considerable, and stratifica- 
tion due to gravity is significant. As a result, the 
medium is dense in its lower heights, but it becomes 
rare and rare as height increases. The perturbed veloc- 
ity amplitude does increase in contrast to the behav- 
ior of perturbed magnetic field. This means that with 




Fig. 4 (Color online) The perturbed velocity mx — z space 
is presented. The panels from top to bottom correspond to 
t = lOr, t = 30r, and t = lOOr, respectively. 



Fig. 5 (Color online) The same as in Figure 
perturbed magnetic field. 



[4] but for the 
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Fig. 6 Time variations of normalized kinetic energy, mag- 
netic energy, and total energy for d = 0.3a are presented 
from top to bottom respectively. 



an increase in height, amplitude of velocity oscillations 
is expanded due to stratification. It is worth to note 
that the density stratification influence on the magnetic 
field is negligible. The divergent configuration of initial 
magnetic field with sheared plasma flow can fasten the 
phase mixing and dissipation of standing Alfven waves 
in Spicules. This is in agreement with the fact that 
spicules have short lifetimes, and are disappeared after 
a few periods. 
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